Macroscopic Quantum Criticality in a Circuit QED 
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Cavity quantum electrodynamic (QED) is studied for two strongly-coupled charge qubits interact- 
ing with a single-mode quantized field, which is provided by a on-chip transmission line resonator. 
We analyze the dressed state structure of this superconducting circuit QED system and the se- 
lection rules of electromagnetic-induced transitions between any two of these dressed states. Its 
macroscopic quantum criticality, in the form of ground state level crossing, is also analyzed, re- 
sulting from competition between the Ising-type inter-qubit coupling and the controllable on-site 
potentials. 

PACS numbers: 42.50.Pq, 05.70.Jk, 85.25.Dq 



Introduction. — In cavity quantum electrodynamics 
(QED) 1], the coupling effects of two atoms Q inside 
a cavity have been theoretically studied to explore exotic 
quantum coherent phenomena e.g., coherent popula- 
tion trapping and dark states. However, the weak dipole- 
dipole interaction between atoms make it difficult to ex- 
perimentally demonstrate these phenomena. 

Recent experiments using on-chip superconducting 
qubits 0, 0, show the strong-coupling effects be- 
tween: i) two superconducting qubits 0- @ : h) a charge 
qubit and a superconducting transmission line resonator 
(TLR) B, @ ; and hi) a flux qubit and an LC oscilla- 
tor llfj . The later two ones have been referred to as 
circuit QED. Moreover, recent experiments |8| demon- 
strated Rabi splitting and AC Stark shift in these QED 
circuits. Combining these approaches we pro- 

pose how to observe the above mentioned quantum co- 
herence effects using two coupled artificial atoms, inside 
a cavity, instead of two natural atoms • 

Here, we present a circuit QED architecture for two 
capacitively-coupled charge qubits, interacting with a 
single-mode quantized field. We not only study the influ- 
ence of strong inter-qubit coupling on quantum coherent 
effects, but also explore its macroscopic quantum crit- 
icality, via level-crossing, for quantum phase transition 
(QPT) 11]. Here, the non-analyticity of the ground state 
is depicted by the fact that the eigenstates are indepen- 
dent of the inter-qubit coupling, while the correspond- 
ing eigenvalues depend on the inter-qubit coupling. Of 
course, a rigorous QPT can only be realized in the ther- 
modynamic limit, but some of the basic features of QPT 
can still be demonstrated in the form of level crossings for 
a system of few qubits 

EJ03 A toy-system for QPT, 
with two qubits, has been experimentally studied using 
NMR 0. 

In contrast with a previous investigation for QPT , 
the two Ising-type-coupled charge qubits in our macro- 
scopic quantum system are also coupled to a TLR. Thus, 



FIG. 1: (color online). Schematic setup of a circuit QED 
with two capacitively-coupled charge qubits, interacting with 
a transmission line resonator. The charge qubits are placed 
between superconducting lines and located at the antinodes 
of a single-mode magnetic field. The couplings between two 
qubits and a single-mode quantized field are realized by the 
magnetic field through the qubit loops. 



for the first time, we incorporate the dressed-state struc- 
ture in the QPTs. Here, we consider a minimal circuit 
QED model to simulate the quantum critical behavior of 
the ground state. 

Circuit QED model with two qubits. — Our circuit QED 
system is schematically shown in Fig. ^ Two iden- 
tical SQUID-based charge qubits are capacitively cou- 
pled to each other with the coupling strength J = 
e 2 C m / (2Cf; — 2C^) through a capacitance C m , where 
Cs is the sum of the capacitances connected to a sin- 
gle Cooper pair box. The level-spacing of each qubit is 
u) a = 2E C (CgVg — 1/2), with capacitance C g , bias volt- 
age V g , and E c = 2e 2 Cs/ (Cf — C 2 J . The two qubits are 
coupled to a single-mode quantized field which is realized 
as a ID TLR with resonant frequency lu — uqtt / (Lylc) . 
Here no and L are the mode number of the resonant mode 
and the length of the TLR; I and c are the inductance and 
capacitance per unit length of the TLR. When the two dc 
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SQUIDs are placed at x n = nL/no, the coupling between 
the qubit and TLR is induced by the quantized magnetic 
field threading the dc SQUID. The model Hamiltonian 
H = Hq + He includes the Ising part 

H Q = ^ a (aP+4V)+J*PaW (1) 
and the Jaynes-Cummings terms 
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H.c. 



(2) 



Hereafter, h = 1, the Pauli matrices a z = 1 1) (T I — 1 1) (I |> 
(7+ = | T) (I I an( l cr - = I 4) (T I are defined by the charge 
eigenstates 1 1) and | J.) denoting 1 or excess Cooper pair 
state, respectively, (a) is the creation (annihilation) 
operator of the resonant mode. The coupling strength 
is g = SEjVfrl<-o/($odVL) with the tunnelling energy 
Ej, the enclosed area S of the dc SQUID, the distance 
d between the dc SQUID and the transmission line, and 
the flux quantum $o = h/2e. 

The symmetry of two qubits enables us to rewrite the 
Hamiltonian H as a function of the total spin opera- 
tors S = + ~aW)/2 and S± = + ct£ 2) , e.g., 

aPai 2) = 2S 2 - 1. The dynamical symmetry SO (3) 
of the Hamiltonian H results in a direct sum decom- 
position of the two-qubit Hilbert space V, spanned by 
{| TT), I 44), I 44), I 44)}; i-e-, V = V<°> © U«. Here, 
V<°> is spanned by the singlet = (| T4) - I 44))/V2, 
while is spanned by the triplet | 44), = (I 44 

) + | 44))/\/2, and | H). and are invariant 

under H. Thus, the total Hamiltonian H can be de- 
composed into a quasi-diagonal matrix with two blocks. 
Obviously, the symmetric couplings of the two qubits to 
the quantized field do not induce transitions between 
the singlet in and any other states in the 

space V^K This is because the collective spin opera- 
tor S± = a± + a± can only change the state vectors 
within an irreducible subspace. Therefore, here, we de- 
note the singlet as a "dark state" 0. This consider- 
ation, based on group representations, automatically pre- 
dicts the coherent population trapping in this supercon- 
ducting macroscopic quantum system, without resorting 
to any dynamical evolution calculations for the natural 
atoms |2j. The population on the singlet \ip~) will be 
trapped to keep its initial value, due to the coherent can- 
celling of the two transitions from both, | fj.) and | J,j), 
to any state. 

Dressed Spectrum. — For quantum criticality, the Ising- 
type Hamiltonian Hq has been extensively studied in 
the thermodynamic limit, i.e., the case with infinite 
qubits To demonstrate this with a few qubits, we 
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consider the eigenvalues of Hq: 
J + uj a , E\y-) = — J and E^^ = J — uj a . We use the 
dimensionless parameter £ = uj a / J to describe the QPT 
character of the ground state. If £ changes from region 



£ < -2 to -2 < £ < 2, and then to £ > 2, the ground 
state of Hq changes from the state | 44) to the state \ip + ) 
(l^ - )) and then to | J.J.). Correspondingly, the "order" 
of the magnetic system changes from ferromagnetic to 
antifcrromagnetic through the maximal entangled state 
\ip + ) {\ip~)) corresponding to the point £ = 0. 

For a weak perturbation, not commuting with 5*2 or 
Hq, the ground state is still determined by the longitu- 
dinal field controlled by the gate voltages. Now, we con- 
sider the effect of the strong interaction on the ground 
state when the Hq system interacts with a single-mode 
quantized field. The above discussions show that the 
quantized field can only mix qubit states, in the invari- 
ant subspace or U' 1 ', to form dressed states, but it 
cannot induce the transitions from to There- 
fore, in the qubit subspace V^ \ the eigenstates \tpn) of 
the Hamiltonian H are the product states of the singlet 
\ip~) and photon number states |n), i.e., \ifn^} = \n,ip~), 

(s) 

which correspond to the eigenvalues En = nuj — J, with 
n = 0, 1, 2,- • •. 

However, in the qubit subspace V^', when we con- 
sider the state mixture induced by the quantized field, 
the Hamiltonian H possesses the following invariant sub- 
spaces : {|0, — 1) = |0,44)}, : {|0,0) = 
|0,^+), |l,-l) = |l,U)},W (n+1) :{|n-l,l) = |n-l>® 
| TT), k 0) = |n) <8 \n + 1, -1) = \n + 1) ® | 44)} 
for n = 1, 2, 3, • • • . Then, the Hamiltonian H can be di- 
agonalizcd in each quasi- diagonal block W^. The 1-D 
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|o,-i) 



block only contains the eigenstate \ipQ 

with eigenvalue E^ = J - uj a . The 2-D block is 
spanned by two eigenstates 

= cos(0/2)|O,O>-Bin(0/2)|l,-l>, 



sin(0/2)|O,O) + cos(0/2)|l, -1) (3) 



with eigenvalues E^ = ±yj J 2 + g 2 . The eigenvalues of 
the 3D blocks W^ n+1 \n ^ 1) can be explicitly solved in 
the resonant case u) = Lu a . In this condition, the eigen- 
values are E n ^ = nuj + J and En — nu> ± N n (g), which 
correspond to the eigenstates 



^ 0) ) = 



(±)\ = 




n-lA) 



2n+ 1 



In + 1,-1) 



- 1,1) ± 



l N n (g) T J 
2N n (g) 



71,0) 



2tt n ± (g) 



|n + l,-l), 



(4) 



where N n (g) = 1/J 2 + (2n + 1) g 2 and to± (g) = 
N 2 (g) T JN n (g). 

Quantum criticality due to level crossings. — The eigen- 
values {E n a ^} with a = 0, ±, s; n = 0, 1, 2, • ■ ■ form a 
complete set for the energy spectrum. Fig. |2 shows a few 
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FIG. 2: (color online) . Rescaled eigenvalues / J} for a = 

0, ±, s; n = 0, 1, 2 of dressed "atom" -photon states versus 
the rescaled longitudinal field £ = w/J, in the resonant case 
uj — uj a and for a parameter ratio p/J = 0.5. 

lowest eigenvalues versus the effective longitudinal field 
£ = w/J, when w a = uj and <?/J = 0.5. However, for very 
weak coupling constant g (e.g., g/J = 0.05), some near- 
neighbor eigenvalues are almost degenerate, for instance, 
E[ +) = uj + J^/l + 3(g/J) 2 « c + J = E[°\ » 
and so on. Comparing with the case without the cou- 
pling to the quantized held, i.e., g = 0, the distribution 
of the spectral structure of the dressed two-qubit system 
becomes very complicated. Usually, the cavity field cou- 
pling g can be very small with respect to J and uj. In 
this case, E^ can be approximated by 

= nuj± [J + S(n)]. (5) 

Eq. Q predicts that level shifts 5(n) = (n + l/2)g 2 /2J 
can be enhanced by the photon number of the quantized 
field. 

To study quantum criticality, we now analyze the 
distribution of the level-crossing points. We plot the 
rescaled 20 lowest eigenvalues En ^/J in Fig. |3{a) for 
n = 0, 1, • • • , 19. It can be seen that all the energy lev- 
els approximately meet together at a certain critical point 
A, which is an intrinsic quantum critical point. This 
is because in the weak-cavity-field coupling limit, i.e., 
g <C J, w, all energy levels E n (n = 0, 1, • • • ) linearly 
depend on n, and thus they are approximately degener- 
ate at a fixed point £o = g 2 / J 2 , independent of n. This 
feature is further demonstrated in Fig.^Jb). We plot the 
eigenvalues of the level crossing points for each pair of 
near-neighbor E n and -E^+i (n = 0, ■ • • , 19) for dif- 
ferent transverse coupling constants g. Fig. |3Jb) shows 
that the crossing points of different levels will slightly 
change along a certain curve for a large coupling con- 
stant g. However, when g gets smaller, the distribution 
of those crossing points gets more compact. In the limit 
of very weak cavity-field-coupling, those crossing points 
converge to the fixed point A with £o = 9 2 / J 2 ■ 



Therefore, when g <C J, uj, the energy spectrum can 
be divided into three parts (shown in Fig. Efb)): (I) 
< £ < £ ; (II) £o < £ < £i; and (III) £ > £ 1; where 
£i = 1 + VI + £o is a degenerate point of E^ and Eq. 
In the regions of (II) and (III), the system has definite 
ground states \G^) = |<4 _) ) and \G^) = \uj { ° ] ), re- 
spectively. It is very exotic that the system does not have 
a ground state in the region (I). This is because E^~' has 
no lower bound in the region (I). Thus, the critical point 
£ = £o is an intrinsic singular point, which is different 
from the generic critical point £i. The existence of this 
intrinsic singular point £o can be verified via the spec- 
trum, generated by the transitions from excited states to 
the ground state. As seen from Fig. |3Ja), the spectrum 
is almost continuous when £ approaches £o from the side 
of £ > £o. The discrete spectrum appears when £ is far- 
away from £o . This feature serves as an experimental way 
to detect the intrinsic singular point in the weak coupling 
limit, as mentioned above. Though we concentrate on the 
resonant case in the above discussions, the results can be 
generalized to the off-resonant case. 

Level transitions and widths of spectral lines. — Above, 
we showed that quantum critical phenomenon can be 
observed from the transitions between different dressed 
states. Now, we further study the optical selection rules 
for the above different dressed energy levels. We con- 
sider that the dressed system interacts with a multi-mode 
bath, then the interaction Hamiltonian H p between the 
dressed system and bath is 

^=E(^ ) 4 1) +5f4 2) )^+H.c. (6) 

fe 

The free Hamiltonian Hi, of the multi-mode bath is 
Hf, = ^ k huJk<i k ak. a\ (&fe) are the creation (annihila- 
tion) operators of the fcth bath mode with the angular 
frequency uJk- The bath can be the quasi- normal modes 
of the TLR or the multi-mode external electromagnetic 
field. 

When g^ k = g k 2 \ H p cannot induce the transitions 
between two different diagonal blocks and of 
the Hamiltonian H; it can only induce transitions be- 
tween different states in the same block. However, when 
g k 7^ g k 2 ^ , the perturbation H, p will break the original 
invariant subspaces and Transitions between 

the different subspaces and are possible. This 
leads to the mixture of the different diagonal blocks of 
the Hamiltonian H discussed above. 

As known in conventional cavity QED, many typical 
strong-coupling phenomena, e.g., the Rabi splitting, ac- 
tually refer to transitions between the different invariant 
subspaces. The corresponding width of the spectral line 
can be determined by the matrix elements ((p„ \H P 
through the Fermi golden rule. We can analyze this prob- 
lem in detail in different quantum critical regions. For 
example, in region (III), |<^o°^) ^ s ^ ne ground state, and 
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FIG. 3: (color online). Quantum criticality described by level- 
crossing points, (a) The rescaled 20 energy levels E„ / J 
(n = 0, 1, ••• , 19) versus £ = u/J with g/J = 0.5. The 
blue line denotes Eq/J = v/I — (g/J) 2 ; red lines show 
E^/J, with n ^ 0. They converge to a critical point j4. 
As a reference, the green line is for E^/J = l-(u/J). (b) 
The same energy levels EnT^/J as in (a) versus the crossing 
points £ = ui/J for near- neighbor levels with different cou- 
plings (g = 0.7J, 0.5J, 0.16J) of qubits to the cavity field. 



\(Pq ') and \cp[ ') are the first two excited states. Tran- 
sitions from \ipo~^) and |<Pg S ') to 1 0, — 1) will lead to a 
double-peak Rabi split. The distance between the center 
of the two peaks is determined by the energy difference 
\J J 2 + g 2 — J between two excited states. For instance, 
if we take J - 4 GHz and g ~ 2 GHz [ljj, this Rabi 
splitting is about 0.47 GHz; however, for a small cav- 
ity field coupling constant, e.g., g ~ 0.2 GHz, this Rabi 
splitting is about 5 MHz. 

We can also evaluate the ratios of the line widths, 
which are determined by the damping rates. For exam- 
ple, the ratio of the damping rates 71 (from |<^q to 
lv?o 0) )) and 72 (from (|^ s) ) to |(^ 0) )) is 



71 (^1) 
72(^2) 



_ G + {uJi)pi{ui) 
M " 2) G-(u; 2 )p 2 (u J2 ) 



(7) 



where pi(u>) is the given spectral density of the bath, 
LJi = w - (2- l)(J + V^ 2 + 9 2 ) (l = 1,2) and G ± (u k ) = 
g^ ± g£ . The above quantitative results can also be 
tested by future experiments. 

Conclusions. — We analyze the dynamical symmetry 
of two strongly-coupled charge qubits, interacting with 
a single-mode quantized field. There exists a coher- 
ent trapped two-qubit singlet state ("dark state"). This 
"dark state" is not affected by any other states when two 
qubits are symmetrically coupled to both the bath and 



the quantized field. However, when the symmetry is bro- 
ken, the "dark state" is no longer "dark" . The transitions 
from this "dark state" to other states become possible. 
By analyzing the Rabi splitting, the level shift in the 
limit of weak-cavity-field coupling, and the dressed state 
structure of the transmission spectrum, we can probe the 
coherent coupling effect between the two-qubit system 
and a single-mode quantized field. 

We also study the influence of the strong inter- "atom" 
coupling and atoms-to-quantized-field couplings on the 
quantum criticality through level crossing. Specially, an 
intrinsic singular point is found in the limit of weak- 
cavity-field coupling. This point is characterized by the 
discreteness of the spectra in some critical regions. Our 
study can be easily generalized to the case of the non- 
resonant interaction between the cavity field and two 
qubits. We hope that our proposal can further moti- 
vate experiments on the circuit QED with two strongly 
coupled qubits. 
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